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We investigate the relation between two-time, multi-spin, correlation and response functions in 
the non-equilibrium critical dynamics of Ising models in d = 1 and d = 2 spatial dimensions. In 
these non-equilibrium situations, the fluctuation-dissipation theorem (FDT) is not satisfied. We find 
FDT 'violations' qualitatively similar to those reported in various glassy materials, but quantitatively 
dependent on the chosen observable, in contrast to the results obtained in infinite-range glass models. 
Nevertheless, all FDT violations can be understood by considering separately the contributions 
from large wavevectors, which are at quasi-equilibrium and obey FDT, and from small wavevectors 
where a generalized FDT holds with a non-trivial fluctuation-dissipation ratio X°° . In d = I, we 
get X°° = i for spin observables, which measure the orientation of domains, while X°° = for 
observables that are sensitive to the domain-wall motion. Numerical simulations in d = 2 reveal a 
unique X°° ~ 0.34 for all observables. Measurement protocols for X°° are discussed in detail. Our 
results suggest that the definition of an effective temperature T c g = T/X°° for large length scales 
is generically possible in non-equilibrium critical dynamics. 

PACS numbers: 05.70.Ln, 75.40.Gb, 75.40.Mg 

Since the analytical solution of the non-equilibrium dynamics of the spherical p-spin model in its 
low temperature phase [l| , many studies have focused on the properties of two-time non-equilibrium 
correlation and response functions, and the relationship between them @, d. In this paper, we 
report on analytical and numerical investigations of several two-time, multi-spin, correlation and 
response functions in the non-equilibrium critical dynamics of Ising models in d = 1 and d = 2 spatial 
dimensions. 

Our work is motivated by the following observations. Multi-point dynamical functions are stan- 
dard objects in equilibrium statistical mechanics which reveal microscopic information related to 
experimentally observable quantities. In non-equilibrated systems, however, the equilibrium relation 
between response and correlation, i.e. the fluctuation-dissipation theorem (FDT), is not satisfied. 
This evident observation became important when it was realized that in the p-spin modelllj] and 
more generally in infinite-range glass models, a generalized FDT can be formulated 0, 0- Q 0| • This 
amounts to the introduction of a fluctuation-dissipation ratio X or, alternatively, of an effective tem- 
perature T e ff = T/X for the slow, non-equilibrated modes of the system Jjjj. The properties of X and 
T e ff have attracted much interest since they suggest that a generalized statistical mechanics can be 
devised to deal with a broad class of non-equilibrium phenomena. 

The generalized FDT is exact for infinite-range glass models only. It is, however, tempting to apply 
the same concepts in other contexts such as glassy systems with finite interaction range, as observed 
experimentally or simulated numerically. A further step is made when those ideas are transferred 
to other physical situations such as domain growth processes 0, @ m non-disordered systems or the 
rheology of soft glassy materials [l0( . Although one does not expect the results for infinite-range 
glass models to apply exactly in all these non-equilibrium situations, it is worthwhile to understand 
analogies and differences, and thus to push these ideas as far as possible. This is the philosophy of our 
paper where non-equilibrium dynamics at criticality is analyzed along the lines described above. Our 
results suggest that the concept of a generalized FDT is indeed useful at criticality, and we describe 
in detail the form it takes as compared to the results obtained in infinite-range glass models. 

The manuscript is organized as follows. The first section below reviews the results obtained for 
correlation and response functions in ferromagnets and delineates the scope of the paper. In Sec. [H] 
the Id Ising model is studied analytically at T c — 0. In Sec. IIIII numerical results for the 2d Ising 
model at T c are presented. A summary and a physical discussion of the results can be found in 
Sec. E| 
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I. FDT AND FERROMAGNETS 
A. Correlation and response functions 



Pure ferromagnets are generally not described as glassy materials, which are loosely defined as 
systems with large relaxation times. However, if a ferromagnet initially prepared at high temperature 
is suddenly quenched to its low temperature ferromagnetic phase, its equilibration time diverges with 
system size llj. This is true also when the quench is performed precisely to the critical point, T = T c . 
In both cases the system remains, in the thermodynamic limit, forever in a non-equilibrated, non- 
stationary state: it exhibits aging. Therefore, one can study physical situations in pure ferromagnets 
that are reminiscent of aging phenomena observed e.g. in spin glasses, polymers or colloids. One 
is then led to ask if the tools used in the glass literature are also useful to describe this type of 
non-equilibrium situation. 

These tools include, in particular, two-time correlation and response functions. Consider two 
physical observables A(t) and B(t). Their connected cross-correlation is defined by 

C(t,t w ) = (A(t)B(t w )) - (A(t))(B(t w )), (1) 

while the conjugate response function is given by 

5{A(t)) 



R(t,t w ) = T 



Sh B {t w ) 



(2) 

h B =0 



Here Kb is the thermodynamically conjugate field to the observable B; for later convenience we 
scale the response by T. Numerically or experimentally it is often more convenient to measure the 
integrated response function, or susceptibility, 

X(*,*w)=y drR(t,T), (3) 

which gives the response to a small constant field hs switched on at the 'waiting time' t w . 

At equilibrium, correlation and response are time-translation invariant and related by the FDT 

R(t,t w ) = ^-C((,(,). (4) 

In that case, (. . .) in Eqs. 111211 stands for the usual ensemble average. If one follows instead the 
dynamics of the system after a sudden quench, the system is out of equilibrium and neither time- 
translation invariance nor FDT are satisfied. Then (. . .) is to be read as an average over initial 
conditions and any stochasticity in the dynamics. In infinite-range glass models, a generalized FDT 
is satisfied in the aging dynamics. The generalization amounts to the introduction of a fluctuation- 
dissipation ratio (FDR), X(t,t w ), through 

d d 
-^-X(t,t w ) = R(t,t w ) = X(t,t w )—C{t,t w ). (5) 
at w at w 

This definition becomes non-trivial because in the limit of long times t and t w the FDR reduces to a 
function of a single variable only, namely the correlation function, 

X{t,t w )^X(C(t,t w )), (6) 

where by a an abuse of notation we retain the same symbol for the FDR and its long-time limit. As 
in equilibrium, response and correlation are then no longer independent quantities, although their 
relationship is now more complex. 

It is now standard to study this generalized FDT via the parametric representation, or 'FD plot', 
of the susceptibility \ as a function of the correlation C 4]. At equilibrium, one has x(i, i w ) = 
C(t,t) — C(t,t w ). Hence, a plot of x(i,t w ) against C(t,t) — C(t,t w ) gives a simple straight line of 
slope 1: this is the equilibrium FD plot. Out of equilibrium, Eq. @ implies 

,C(M) 

x(Mw)= / dxX(x). (7) 
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Therefore, when Eq. JJ]J holds, the FDR can be obtained directly from the slope of the FD plot, 
which is X(C). Otherwise, from Eq. JSJ a plot of x(t, t w ) against C(t,t) — C(t,t w ), with t fixed and 
t w the curve parameter, will still have slope X(t, t w ). Since the amplitude of correlation and response 
functions can diverge or converge to zero for t — > oo (see below) it can be useful to use normalized 
quantities, plotting x(t,t w ) = x(t,t w )/C(t,t) versus 1 - C{t,t w ) where C(t,t w ) — C(t,t w )/C(t,t). 
Since the normalization factors are independent of t w , the slope of the FD plot is then still given by X. 
The normalization issue is less important when presenting numerical data, which are by construction 
obtained in a restricted time window where the amplitudes of the dynamical quantities typically 
change only slowly. 

Appealingly, the FDR can also be interpreted as defining an effective temperature, T c g(i,t w ) = 
T / X(t,t w ), replacing the equilibrium temperature T by an equivalent quantity out of equilibrium. 
Moreover, it is a general result that for the case A = B, where one considers the autocorrelation 
of A and the associated response, T a e is the temperature measured by a thermometer coupled to 
the observable A at the appropriate timescale .6]. As a direct corollary, this effective temperature 
then satisfies the zeroth law of thermodynamics. Clearly, however, the introduction of an effective 
temperature is of thermodynamic interest only if T e s is actually independent of the observables A 
and B under consideration. This is indeed true for infinite-range glass models @, implying that 
although the system is out of equilibrium it can still be described in thermodynamic terms, at the 
moderate cost of introducing one extra parameter, namely the effective temperature Beyond 
infinite-range glass models, the observable-dependence of the effective temperature remains largely 
an open question but has been discussed in detail in the context of trap models [13J and in a realistic 
numerical model of a supercooled liquid . 



B. Ferromagnets at low temperature 

For ferromagnets, two-time dynamical functions have been studied both for a quench to the low 
temperature phase and to criticality, with most work to date focused on the first situation. 

In the low temperature phase, the evolution of the system consists in the growth of ordered domains, 
with a typical domain size £{t). Two-time quantities that have been thoroughly studied are the spin 
autocorrelation function, C s (t,t w ) = (si(t)si(t w )) , where Si(t) is the value of the spin at site i at 
time t, and the conjugate response function, R s (t,t w ) = — {d/dt w )x s (t,t w ) = 5(si(t)) /Shi(t w ), where 
hi is the magnetic field at site i. In this case, the connected and disconnected correlations coincide 
since (s Jt) ) = at all stages of the coarsening process. From the analytical solution of solvable 
models [H 01 111 111 El and the simulation of more realistic situations 011 E3, El E2, the 
behavior of these two quantities is now well known. For small time differences, At = t — t w <C t w , 
time-translation invariant behavior is observed, C s (t,t w ) « C s (At), Xs(Mw) ~ Xs(At), and the FDT 
Xs(Mw) ~ C s (t,t) — C s (t,t w ) is obeyed. This first regime is due to thermal fluctuations in the bulk 
of the domains, which are essentially equilibrium fluctuations. For larger time separations, At 3> t w , 
the fluctuations of the interfaces dominate the dynamical behavior. The reasonable hypothesis that 
coarsening is a self-similar process, in the sense that for large times all dynamical functions depend on 
time only through the typical domain size £, implies the scaling form C s (t, t w ) ~ C s (£(t)/£(t w )) for the 
correlation function. The contribution of the interfaces to the response function can be estimated 
as a sum over all wavevectors, Xs(t, ^w) ~ Jy^u \ d d k fc" 2 £ _1 (t). This expression results from the fact 
that the response at time i w is dominated by large wavevectors, k£(t w ) ^> 1, each wavevector k giving 
a contribution of the order of k~ 2 . The factor £~ 1 (t) represents the density of domain walls, and a is 
a UV cutoff given by the lattice spacing. This reasoning implies for the long-time contribution to the 
susceptibility the scaling form Xs{t 7 t w ) w f(£(t w ))xs{£(t)/£(t VJ )), where the function f(x) depends 
on the dimensionality d of space and is given by f(x) = 1/x for d > 2, f(x) = (\nx)/x for d = 2 
and f(x) — const for d = 1. This scaling function has recently been revisited in Refs. jHI^, with 
particular attention to the case d = 2 [3 H3 • 

From the above arguments, and for d > 1, the parametric plot of Xs(t, t w ) versus C s (t, t w ) — C s (t, t w ) 
consists of an initial equilibrium part followed by an essentially horizontal section. In the latter 
the correlation function decays due to interface motion, while the response function hardly changes 
because any contribution from the interfaces is suppressed by the f{£(t w )) prefactor. If a limiting 
FDR, A s °°, is defined through 

A s °° = lim lim X s (t,t w ), (8) 
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then it follows that X^° — for d > 1 in coarsening processes. For d = 1, on the other hand, both Xs 
and C s are scaling functions of £(i)/£(t w ) and the parametric plot assumes no simple shape, implying 
that X,? could be any finite number. This is confirmed by the analytical solution of the dynamics 
of the Ising chain at T = which shows that A s °° = \ j^l H3|. The factor i was first derived in 
Ref. 0. 

C. Ferromagnets at the critical point 

The non-trivial value of X£° for the Ising chain was interpreted u sing the fact that in d = 1, the 
ordering temperature T = coincides with the critical point, Tn= |2g. It was then suggested that 
a non-trivial X£° could be a generic feature of critical points |29j . This is physically reasonable, since 
the whole argument for = in coarsening processes relies on the separation between bulk and 
interfaces; this is no longer valid at the critical point where the bulk has the well-known self-similar 
structure of ferromagnets at criticality. 

From analytical and numerical studies, the behavior of two-time single-spin dynamical quanti- 
ties is again well understood, as briefly reviewed in Ref. |3f| . Physically, the non-equilibrium 
dynamics following a quench to the critical point consists in the growth of the dynamical corre- 
lation length, « t 1 ^, where z is the dynamical critical exponent j^l). Critical fluctuations 
of large wavevectors, k£(t) 3> 1, are almost equilibrated, while those with small wavevectors, 
k£(t) <C 1, still retain their non-equilibrium initial condition. This separation leads to the seal- 
ing forms C s (t,t w ) « Ai~ 2 ^ z C s (^)/^w)) and x a (t,t w ) » Ai- 2 ^ Xs (£(t)/£(t w )), where p and 
v are the standard critical exponents. This can be interpreted as follows. For short time differences 
Ai <C t w , equilibrated fluctuations with large k dominate and dynamical functions assume their stan- 
dard equilibrium power-law decay. The dynamics at large time separation At 3> t w , on the other 
hand, is dominated by the growth of the dynamic correlation length and leads to the £,(t)/£,(t w ) scal- 
ing. This in turn implies that, beyond the initial equilibrium part, the FD plot will again assume 
a non-trivial shape, as in the Ising chain. The striking similarity of these results with the aging 
dynamics of finite dimensional spin glasses was noted in Refs. (2^, ■ 

The reasoning above confirms that at criticality X£° can take any finite value, in contrast with the 
A s °° = obtained in the low temperature phase. It was further argued that X^° should be a new 
universal quantity at criticality 29] . As such, it can be computed using standard renormalization 
group procedures, and this program has recently been started for various models [33|, |]3 E3 The 
value of X s °° is known exactly for the Ising chain [2^, l2ll l28j . where A s °° = |, for the spherical 
ferromagnetic model [2j|, where X^° — | for d > 4 and A s °° — | for d = 3, and for the Gaussian 
model 0| where X£° = | . An estimate is known for model A at second order in 4 — d , to first 
order in y/4 — d for the diluted Ising model [34| , and to first order in 4 — d in model C j3^ . 

D. Motivations for this work 

This short review of known results in the non-equilibrium dynamics of pure ferromagnets shows 
that much research has been done on the subject. So, why another paper? 

First of all, the relevance of the notion of an effective temperature at criticality can be questioned 
because the FD plots for the spin dynamic functions do not assume a simple linear shape with a 
well-defined slope, as happens in the low temperature phase. This is related to the fact that at low 
temperatures the decay of correlation functions occurs on two well-separated time scales. Each has 
its own associated effective temperature, a fact reminiscent of the physics of structural glasses. At 
criticality, on the other hand, one has a continuum of time scales associated with different wavevectors, 
t(k) ~ k~ z . Moreover, for finite k the equilibration time is finite, meaning that the number of modes 
that are still out of equilibrium decreases as time increases. This suggests that an effective temperature 
could be relevant only when considering the k — * limit, a point which our analysis will clarify. 

Second, we mentioned the important issue of the observable-dependence of a generalized FDT. 
This issue remains completely open since the studies cited above focused exclusively on a single 
FD relation, for the spin autocorrelation and associated response. In order to get a more complete 
theoretical understanding it is crucial to understand if other observables give the same results, and if 
not, how they are related. 
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A third motivation for the study of higher order correlation functions comes from the observation 
that the dynamics of coarsening models is dominated by the motion of topological defects. For 
Ising models, these are domain walls, the local density of which is given by the 'defect' observable 
Si(t)sj(t), where are nearest neighbors. Defect dynamical functions have recently been studied 
in the context of kinetically constrained Ising models j3|J , and the FD relations that arose showed 
interesting and unexpected features. 

For an Ising model, there are at least four 'natural' FD relations, involving respectively the spin 
autocorrelation, C s (t,t w ), the magnetization density m{t) correlation, C m (t,t w ) = (m(t)m(t w )) , the 
defect autocorrelation, C d (t,t w ) = (si(t)sj(t)si(t w )sj(t w )) — (si(t)sj(t))(si(t v ,)sj(t w )) with near- 
est neighbors, and the energy density e(t) correlation, C c (t,t w ) = (e(£)e(i w )) — (e(t))(e(f w ))- Note 
again that connected and disconnected correlation functions coincide for the magnetization; this is 
not the case for C d (t,t w ) and C c (t,t w ). In the Id case, we will also investigate two-time functions 
which smoothly interpolate between incoherent, local functions (spin or defect) and coherent, global 
ones (magnetization, energy), and discuss the case of correlation functions of higher order. In the 2d 
case, we will stick to the four quantities listed above. 



II. THE Id ISING MODEL 



In this section we study the non-equilibrium dynamics in the Glauber-Ising chain with Hamiltonian 



H = ~y^Sj flj+i, (9) 

i 

where the Si (i = 1,...N) are N Ising spins subject to periodic boundary conditions. Glauber 
dynamics consists in each spin Sj flipping with rate ^[1 — ^Si{si~\ + Sj+i)], where 7 = tanh(2/T). 
We focus on the evolution of arbitrary two-time spin and defect correlation and response functions 
in the thermodynamic limit N — > 00, after a quench from equilibrium at T = 00 to T — > 0. As 
explained above, although a variety of aspects of the associated coarsening dynamics have already 
been studied l37| , results on the non-equilibrium FDT violation so far are restricted to the spin 
autocorrelation and response functions [2(jl |2?| . In Sec. Ill Al we introduce the more general class of 
spin and defect observables we investigate. We briefly present the main result of our method to derive 
multi-spin two-time correlation and response functions, as developed in Ref. [3S| ]. and summarize the 
approach used to extract from this the quantities of interest here. Our results for spin dynamical 
functions arc then given in Sec. Ill Bl — a preliminary account of which has appeared in the conference 
proceedings — and for defect functions in Sec. Ill CI In Sec. IIIDl we discuss the physical aspects 
of our results for the Id Ising model. 



A. General strategy for the calculations 



We consider the following spin and defect observables O s and Od 



(10) 



In both cases the ej arc quenched random variables with zero mean [ej] = and translation invariant 
covariances q^j = [ciCj]; here [•] denotes the average over the distribution of e. Without loss of 
generality we set qo = 1. We define the corresponding connected two-time correlation functions 

C(t,U) = ^[(O s (t)O s (t w ))} and C(t,t w ) = ^[(O d {t)O d (t Vf ))-(O d (t))(O d (t w ))], (11) 



for spins and defects, respectively, and the responses 



T 



S(O s (t)) 



Sh s (t v 



h.=0 



T 

and R(t, t w ) = — 



Sh d (t w ) 



(12) 



h d =0 



where h s and h d are thermodynamically conjugate to O s and Od, respectively. All functions are scaled 
by N to get quantities of order unity. It is easy to show that, in the thermodynamic limit N — > 00, 
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Eqs. d), C3 become 38] 

C(t,t w ) = y^ j q n Cn{t,Ut) and R{t,t w ) =^q n R n (t,t w ). (13) 

n n 

Here we have used translational invariance (which holds for our quench from an equilibrium state) to 
define the distance-dependent correlation functions 

Cj-i(t,t w ) = (si(t)sj(t w )) (spins), ^ 

Cj-i(t,t w ) = (s i (t)si + i(t)s i (t w )sj + i(t w )) - (si(t)s i+ i(t))(sj(t w )sj + i(t vr )) (defects), 

and associated responses 

R j _ i (t,t w )=T S{Si{t)) '-^ — ' - ,^_,%W S .«W) 



8hj(t w ) 



(spins) and Rj-i(t,t w ) = T 
hi=o on jJ+i\ t «) 



(defects). (15) 



1=0 



As usual, hj and hjj+i are conjugate to Sj and SjSj+i, respectively. Translation invariance also 
shows that in the thermodynamic limit N — > oo, the expressions (|11|) . (|12fl are self-averaging, i.e. 
independent of the particular realization of the disorder variables 

Analysis of the non-equilibrium FDR for the observables O s , Od thus requires knowledge of all spin 
and defect correlation and response functions (|14fl . (|15fl . We have tackled this problem in [38| where 
we give closed, exact solutions for generic two-time multi-spin correlation and response functions in 
the Glauber-Ising chain after a quench from an arbitrary equilibrium state at temperature 7] > 
to any T > 0. The approach is based on the closed hierarchy of differential equations for the spin- 
correlations (Sjj (t)si 2 (t) - ■ ■ Si k (t)), which we solved for arbitrary initial correlations (s^ (0) • • • Si k (0)). 
The key result reads [38| 

LU i 

(•*,(*)-■*.(*)> =E E (- 1 ) OT n H Max,-Max- 1 )(2t)*!L-:t,..M,))W ^ 
;=0 ireV(l.k) A=l 

*i fc) (*)= E ( E (-ir]]e-\-J, w (7t)]W(0)-"«A(0)) (17) 

ji<-Jk \7reS(k) A=l / 

Here 7r denotes permutations and (— 1)^ their sign; S(k) is the set of all permutations of {1, 2, . . . k} 
while V(l,k) is the set of permutations corresponding to choosing I ordered pairs from the numbers 
1, 2, . . . k and keeping the remaining k — 21 numbers in ascending order. Explicit expressions for the 
functions and H n (x) are giv en in |3_8] (for N — > oo the I„(x) are modified Bessel functions, see 

Appendix A). We also show in |38j that the evolution of two-time, multi-spin correlation and response 
functions is governed by an identical hierarchy of differential equations, so that these quantities can 
be obtained from ifTBl) . ltT7|) if we substitute the corresponding equal-time initial conditions in (|T7jl . 
The latter are just equal-time correlations — or can be expressed in terms of these for the response 
functions — which we know already. For a quench from an equilibrium state this leads to explicit 
results for the two-time multi-spin functions. As simple examples we state in 38] the spin and defect 
functions l|14|) , CH for the quench from 7] — oo to T — > considered here. For spins one finds 

C„(Mw) = e-( t+t ^{l n (< + < w ) + J At Ut + tw-r) Po+Ii](t)}, (18) 

o 

t 

Xn(t,t w ) = ie"* J dr e- T l n (t-r) po + 2Ii+I 2 ](2r), (19) 



and for defects 

C„(t,t w ) - I e -(' +t -) [!„_!-!„+!] (t + t w ) f dre- T [I n _i-I n+1 ](r) (20) 



2 

t-t 



+ e" 2t 



- *w) + 2I„ + In+iK* + tw) ~ e- 2 *" [(I n _! + I„)(I„ + I n+ i)]{t + t w )}, 

Xn(Mw) = e- 2 *{2^, [I(, + Il](2i)-In(*-*w) [I„-l + 2I n + I n+1 ](i + i w )}. (21) 
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Here and below the short-hand [ . . . ] (x) is used to indicate that all functions enclosed in the square 
brackets have the same argument x; S ny o is the standard Kronecker delta. The expressions 119fl . I(21|) 
for the susceptibilities are more convenient than those for the responses R n (t, f w ) — ~(d/dt w )xn(t, f w ) 
and so we mostly base the following discussion on them. We note that while equations l|18Jl , (|19|l have 
already been given in various forms, e.g. [2(J, we are not aware of any equivalent in the literature of 

Eqs. ((18(1 - 1(21(1 will form the basis for our analysis of the FDR in the Id Ising chain in Sees. Ill Bl 
and III Cl For now we return to the observables O s , Od and in particular the choice of the field 
covariances q n . According to 1(13(1 we obtain spin and defect autocorrelation and response functions 
by choosing uncorrelated random fields i.e. q n = S n fi. We abbreviate the notation in this case to 
that used in the introduction and write C s (t, t w ) for spin and Cd(i, t w ) for defect autocorrelations and 
similarly Xs(t,t w ), Xd(t,t w ) for susceptibilities. Uniform covariances q n = 1, on the other hand, yield 
full summations over all cross-correlation and response functions in lj*H?)l . So O s and Od produce just 
the magnetization and energy, respectively; we thus use the obvious short-hands C m (t,t w ), x m (t,£ w ), 
and C e (t,t w ), \e(t,t w ) for this case. It will turn out that the local {q n — 6 n fl) and global (q n 
I FDT relations for spin and defect observables are very different. Therefore we also investigate 
intermediate choices of q n that interpolate between these two extremes. Two classes of covariances 
can be distinguished. We may interpolate between q n — S n fi and q n — 1 by a family of covariances 
that satisfies J2 n \ln\ < oo for any non-uniform choice of q n ; we call the corresponding fields e, short- 
range correlated. Alternatively, we can interpolate such that ^2 n \q n \ = oo as long as the fields are 
not completely uncorrelated; we refer to such fields as infinite-range correlated. In cither case the 
analysis of the FDR for the correlation and response functions of the associated observable requires us 
to evaluate the infinite sums in ((13(1 . This can be done conveniently in terms of the Fourier transforms 
q(k) = T{q n }, C(fc; t, t w ) = :F{C„(£, t w )} and x(k; t, t w ) = T{x n {t, t w )} where 

TT 

Hfn}=Y,^ e ~ mk and ^H/(*)} = / ^/We 1 ""- (22) 

— TT 

In Appendix IE! we state the Fourier transforms of l|18l) -- (|21[l . in terms of which (|13|1 becomes 

7T 7T 

/dk f dk 

— q(k)C(k;t,t w ) and x(t,t w ) = / — q{k) X (k;t,t w ). (23) 

— 7T —TT 

An explicit example of a family of short-range correlated fields, parameterized by a > 0, is given by 
the Lorentzian covariances 

2 

9l,„ = ^- — o qh(k)= a * cosha(7r- |fc|). (24) 

a 1 + n z smh an 

The transform l|24|l can be found in any table of Fourier transforms. Eq. (|24|1 indeed defines short- 
range correlated fields: since <7L,n > the criterion becomes N c = ^ n qL.n < oo which is satisfied 
since N c = 9l(0) = a7rcotha7r. By varying a we can also smoothly tune our observables between 
local (qh,n — * #n,o as a — * 0) and global (qh,n —> 1 for a — > oo) ones. We denote the corresponding 
correlations and susceptibilities by Cl(£, t w ) and Xh(t,t w ). An example of covariances that yield 
infinite-range correlated fields is 

where < a < 1 and T(x) is the Gamma function It is clear from (|25J) that qp^ n is even in n 
and qpo = 1. We show in Appendix ICl that for a —> 1 we get qp, n — S ni o while a->0 gives qp :7l = 1. 
We also prove there that qp, n decreases monotonically as |n| increases, decaying asymptotically as a 
power-law qp, n ~ I 71 ! -0 ! and that indeed J-~ 1 {qp(k)} — qp in . The reverse transform ^{gp,™} does 
not converge in the usual sense, but this is not necessary for equivalence of 1)13(1 and (|23|l . So 1)25(1 
again allows us to interpolate smoothly between local and global observables, but in such a way that 
^2 n Qp.n = oo for any a G [0, 1[. The correlations and susceptibilities for the observables defined by 
fields ej with the power-law covariances <(25|l are denoted by Cp(t, f w ) and %p(f, t w ) below. 
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B. Spin observables 

1. Random field: Incoherent functions 



The FDT violation for the spin autocorrelation and response function has already been studied in 
detail in pfj. W\ . In particular it was shown that the FD plot approaches a non-trivial limit curve 
in the aging regime, with X^° — \ . We can easily recover the existing results for C s (t,t w ), Xs(Mw) 
from our exact solutions JTS), i|19f) by setting n = 0. It is useful to focus on the aging limit. Formally 
this is an asymptotic expansion in the limit t, t w — > oo with e < t w /t < 1 — S fixed and e, S > 0, to 
ensure that t, t w and At all diverge and are of the same order. In this limit the asymptotic expansion 
13(1 for modified Bessel functions yields immediately 



2 / 2t w 

C s (t,t w ) ~ -arcsinW w , (26) 

7T V t + t w 

Xs(t,t w ) ~ — arccosi/^. (27) 

7T V t 

Here and below, the '~' sign denotes results which are asymptotically exact in the aging limit. The 
limit FD-plot corresponding to ((26(1 , 127(1 is contained in Fig. [21 below and the associated FDR is a 
function of the time ratio t w /t only, 

^ s (Mw)~^. (28) 

It shows a continuous crossover from X s (t, i w ) = 1 for At <C t w to X s (t, t w ) — X^° — | for At 3> 
We note that the aging expansion of the spin correlations and susceptibilities 1(18(1 , l(19(l is dominated 
by the leading term of the asymptotic series l(A3(l for the modified Bessel functions, which is inde- 
pendent of the order n. Therefore 1(26(1 . 1(27(1 in fact apply to all finite-distance spin cross-correlations 
and susceptibilities C n (t,t w ), Xn(t,t w ). Consequently, the latter produce the same limiting FD plot 
and FDR Jl| as for n = 0. 



2. Uniform field: Coherent functions 

As described, the uniform field effectively allows us to study the FDR for the magnetization. 
The corresponding correlation and susceptibility are most conveniently obtained from the Fourier 
transforms ((Bl(l . i(B2(l by setting k = 0; the time integrals appearing in C(k;t,t w ), x(k',t,t w ) can 
then be solved. One finds 

C m (Mw) = e- 2t -{I (2i w ) + 4t w [I + I 1 ](2i w )}, (29) 
# m (Mw) = ^e- 2t -[I + 2I 1 +I 2 ](2t w ). (30) 

We have given the response R — —dx/dt w here rather than the susceptibility because it has a simpler 
form. Note that both the correlation and response function ((29(1 . 13U|) are independent of t. This 
can be understood from the fact that for T = the magnetization m = {l/N)J2i s i performs a 
random walk with step size ±2/iV and a time-dependent rate \Nc, where c is the concentration of 
domain walls. The latter can be obtained explicitly from 1(18(1 by setting t = t w and n = 1, since 
Ci(t w ,t w ) = ( Si (t w )s i+1 (t w )) = 1 - 2c(t w ). This gives 

c(tw) = ie- 2t »[I + I 1 ](2i w ). (31) 

Now, since the random walk of m is unbiased, C m (t,t w ) = N (m(t)m(t w )} — C m (t w ,t w ) follows 
immediately. Also, C m (t w ,t w ) will grow with rate 2ND(t w ), where D = (2/jV) 2 §iVc = 2c/ N is the 
diffusion constant. One should thus have dC m (t w ,t w )/dt w = 4c(i w ), and from (|29|l . (fTTTf) one verifies 
that this is indeed the case. Similar arguments apply to the response R m (t,t w ). Brief application of 
a field at t w biases the domain wall motion and hence the random walk of m; thereafter the random 
walk is again unbiased and so the response must be ^-independent. The momentary bias in the 
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random walk rates contains the domain wall concentration c(£ w ) as an overall factor, and consistent 
with this expectation one gets R m (t,t w ) ~ 2c(£ w ) asymptotically. 

Since both C m (t,t w ) and R m (t,t w ) are functions of t w only this also applies to the FDR 

[I + 2I 1 +I 2 ](2t w ) 
XM - 4Po + I 1 ](2t w ) ' (32) 

which crosses over from the initial value X m (0) = \ to h on an 0(1) time scale. So, apart from a 
transient after the quench, we measure X m (t w ) — i for all t > i w 3> 1; in particular the limiting value 
Am(i w ) ~ -^"m = ^"s^ = 1/2 in the aging regime is identical to that for the incoherent spin observables. 
Note that there is no quasi-equilibrium regime with X m = 1 for At <C f w . The corresponding FD 
plot converges to a straight line of slope | (see Fig. below). 



3. Short-range correlated field 

Next we investigate the effect of short-range correlations in the random fields e, on the FDR. The 
correlations and susceptibilities of the corresponding observables may be obtained either from a real- 
space summation (1131) or an integration in the Fourier representation (|23|1 . Using the latter, we note 
first that the short-range criterion ^ \q n \ < oo for the covariances implies that q(k) = T{qn\ = 
J2 n 1 ne ~ mk i s a continuous function. The Fourier transforms C(k;t,t w ), x(k;t,t w ), on the other 
hand, satisfy 

^2^1^2n~S(k) and -> 2^(Jfe) (33) 

C s (t,t w ) Xs(t,t w ) 

in the aging limit, where <5(-) is a 27r-periodic version of the ordinary Dirac delta. The normalizations 
of the right-hand sides of (|33|l are clear since C s (t,t w ), for instance, is given by the Fourier integral 
f2"2"|l over C(k;t,t w ) for n = 0. And since C s (t,t w ) and Xs(Mw) are 0(1) functions of t w /t in the 
aging limit, while C(k;t,t w ), x(fc;i,i w ) vanish in the same limit for any k that is not a multiple of 
27r, Eq. i|33|) follows. This in turn implies that for any short-range correlated field 

C(Mw) ~ N c C s (t,t w ) and x(*.*w) ~ N cXa (t,t w ) (34) 

provided that N c — q(0) — J^nln, which estimates the number of lattice sizes over which field 
correlations extend, is nonzero. So in the aging limit the correlations and susceptibilities are ultimately 
just proportional to Q26[l. (|27|l and hence yield the same FDR and FD plot as the local spin functions. 
This statement can equivalently be made in real space, based on convergence of the series q n and the 
fact that all finite-distance cross-correlation and response functions behave asymptotically as l)26[l. 

At finite times, however, we find a crossover between two dynamical regimes. A scaling analysis 
shows that the peaks in C{k\ t, t w ) and x(fc; t, t w ) at k = have widths t~ x / 2 and Ai -1 / 2 , respectively. 
Correspondingly, we have growing length scales in real space. These are ic ~ t x l 2 for correlations, 
corresponding to the typical domain size, but i x « At 1 ' 2 for the response which reflects the fact that 
perturbations spread diffusively. When icJ-x one has an effectively local observable and we 

are in the asymptotic regime (13411 . If, however, £cJ-x ^ ^ci the fields e, are correlated over distances 
much longer than the dynamical length scales, giving an effectively uniform field. One thus expects 
to get an FD plot similar to that obtained for the magnetization. The illustration of the crossover in 
Fig-E obtained by numerical integration of l|23ll . shows that this is indeed the case. 



4- Infinite-range correlated field 



For infinite-range correlated fields one cannot use simple scaling arguments since the correlations 
and susceptibilities contain contributions from all length scales. For the power-law covariances l|25(l 
introduced above, this is reflected in the singularity of qp(k) at k = 0. Therefore we have to analyze 
the full expressions for Cp(t,t w ) and xp(^w) that follow from (|2~3l after substitution of gp(fc), 
C(k; t, i w ) and t, t w ). Fortunately, for the particular choice of qp(k) the results may be expressed 
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FIG. 1: Normalized FD plot (left) and the corresponding FDR versus t w /t (right) for a random field with 
Lorentzian covariances (1241 with a = 10 3 , correlated over N c = an coth an « an sites. In the FD plot, t is fixed 
for each curve and varies over the ranee 10°, 10 1 , . . . 10 10 (bottom to top). The lines first converge towards 
the straight line with slope | corresponding to a coherent observable (the magnetization) but eventually, for 
t > 10 , cross over to the limit plot for uncorrelated fields. This behaviour is also reflected in the evolution 
of Xi,(t,t w ). There, however, we have the freedom to fix either t or t w . The plot shows the case of fixed t w , 
which is more convenient for comparison with simulations, for i w = 10°, . . . 10 10 (bottom to top). 

in terms of single integrals of the form 

C P (t, t w ) = e- 2{t+t ^ jiFi (|, *±S; 2(t+t w )) +Jdr e T iFi (§, i±a ; 2(t+t w -r)) [I + Ii](r)| , (35) 
t 

Xp(Mw) = ^ e " 2 '/ dT i F i (h ^ 2 ( < - T )) Po + 2Ii + I 2 ](2r), (36) 

tw 

where i,Fi(a, 7; z) is the confluent hypergeometric function [4C|. Equations. (|35|) . I|36|l are exact and 
can be used to study the FD plots and the FDR numerically. However, in the aging limit asymptotic 
expansions may be substituted for the non-elementary functions and significant simplifications are 
possible. One finds 

Cp(Mw) ~ ^-r(i±s) (t + t w ) i ^B(i,l-f; i f--), (37) 

XP(*,*w) ~ ^ 1 r(l±^) ^ [B(i,l-f)-B(i,l-f;^)] , (38) 

where R(p,q;x) is the incomplete Beta function B(p, q;x) — duu p ~ 1 (l — u) 9_1 and B(p,q) = 
B(p, 5; 1) is the complete one |40j. In the random field limit, a — > 1, we recover the expansions (12611 . 
Ij27(l for the incoherent functions since B(^, ^; x) = 2 arcsin-^/x, whereas the uniform field limit a — > 
can be shown to coincide, using B(i, l;x) = 2y/x, with the asymptotic expansions of (|29|l . (|30|l for 
the coherent functions. So the power-law covariances l|25[l indeed allow us to interpolate between the 
coherent and incoherent observables. For intermediate exponents < a < 1 the fluctuations in the 
observable O s grow as t^~ 2 and the two-time correlation (|37|l has a plateau at a corresponding 
value for At -C t w ; for At » t w it decays as tw~ a ^ 2 (t w /At) Q / 2 . For the susceptibility we deduce 
from QgJ a A*( 1_a )/ 2 (At/*w) 1/2 growth for At < t w that crosses over to At^"'/ 2 for At > t w . 
Fig. shows exact FD limit plots that follow from lpTf)l. (j^BJ. The associated FDR may be obtained 
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from jS gfl as 

In principle one should first differentiate ffify . (fTTF>f) to obtain Rp(t, t w ) and (d/dt w )Cp(t, i w ) and then 
perform the aging expansion, but this turns out to give the same result. Equation (|39|l is a function 
of t w /t only and interpolates between the FDR (|28|l for the local spin observables (a — > 1) and the 
constant = | for the magnetization (a — * 0). Plots of Xp(t,t w ) for various powers a are also 
shown in Fig.El It is remarkable that the FDR again crosses over from Xp(t, t w ) = 1 for At -C i w to 
= I f° r &t^>t w , independently of the power-law exponent a. 




FIG. 2: Normalized FD-plots (left) and the corresponding FDR versus t„/t (right) in the aging limit t,t„ — > 
cxd. The bottom curves in both plots are for the magnetization (Sec. Ill B 20 and coincide with the uniform 
field limit a — > of power-law field covariances. The intermediate curves are for power-laws (Sec. Ill B~4t with 
exponents a = 0.2, 0.4, 0.6, 0.8 (bottom to top). The top curves represent the random field limit a — * 1 for 
power-laws, and apply also to any short-range correlated field (Sec. Ill B 31 with N c / or, in the extreme 
case, the incoherent functions of Sec. Ill B II 



5. Harmonically correlated fields and X 



The explicit examples given in Sec. Ill B 1III B 41 suggest that X°° = ^ is a generic feature for the 
spin-observables O s defined in Hl()(l . To show that this is indeed true, we start from the fact that for a 
general observable O s the correlation and susceptibility — and hence (d/dt w )C(t,t w ) and R(t,t w ) — 
may be written in the form l|23|l . By introducing a generalized FDT for the Fourier modes 

R(k; t, t w ) = X(k; t, t w )^-C(k; t, t w ), (40) 

we may express R(k;t,t w ) via H4()[) and thereby obtain the following representation for the FDR 
X(t,t w ) associated with a generic spin-observable O s : 

J §X(k;t,t w )q(k)^C(k;t,t w ) 
X(t, t w ) = . (41) 

/ §Q(k)m^C(k;t,t w ) 

— TV 

This means that X(t, t w ) may be considered as the average of X(k; t, t w ) over the normalized distri- 
bution of q(k)(d/dt w )C(k;t,t v/ ) on k £ [— ir, ir]. The FDR for Fourier modes follows from (|40|l and 



12 



the expressions I|B1|) . (|B2|l for C(k;t,t w ), x{k;t,t w ) as 



, x I + 2l! +I 2 2t w 

X(k;t w ) = p — T , (42) 

4[I + Ii](2t w ) - 2(1 - cosfc) |e 2t » cosk + J w dre( 2i .- T >™ s ' ! [Io + Ii](r) J 

and is a function of i w and A: only. For fc = 0, l)42|) reduces to the FDR for the magnetization X m (t w ) 
ijS^ and hence X(0; t w ) w i for t w 3> 1. A scaling analysis of shows that for \k\ <C 7r and t w ^> 1 
we get X(k;t w ) « X(k 2 t w j with X(k 2 t w ) w | for £; 2 i w < 1 and X(k 2 t w ) w 1 when fc 2 i w > 1. So 
H42(l reflects the successive equilibration of increasing length scales. 

Now we can return to the FDR l|41l) for the observable O s . For the magnetization — being the 
coherent observable — we have q(k) = 2n5(k) and (|41|l reduces to the trivial identity X m (t w ) — 
X(0; t w ). In physical terms, by selecting the coherent observable we only measure the FDR associated 
with the infinite length scale. For other spin-observables, being characterized by the function q(k), the 
FDR X(t, t w ) contains contributions from all length scales. For the long-time limit X°° , however, the 
situation simplifies because (d/dt w )C(k;t,t w ) develops an infinitely sharp peak at k — as t — ► oo. 
This can be verified by a scaling analysis of (|B1I) . For sufficiently well-behaved functions q(k), the 
normalized version of the distribution q(k)(d/dt w )C(k;t,t w ) thus becomes a realization of 8(k) and 
we get X(t, t w ) — > X(0; t w ) as t — > oo. Taking the limit t w — > oo then shows that X°° = i, as claimed. 
So for a generic spin observable, X°° again just gives the FDR associated with the infinite length 
scale. The only exception occurs when this contribution is explicitly suppressed. An example of the 
latter case would be harmonically correlated fields, q n — cosnp with < p < tt: for such observables 
X(t,U) = X(p;t w ) and hence X°° = 1. 



C. Defect observables 



1. Random field: Incoherent functions 

The defect observable Od given in (jl()|l with random, uncorrelated fields €i allows us to study the 
FDT violation for local defect correlations and susceptibilities. These follow from 1)21) [) . I|21l) by setting 
n = 0, giving 

C d (t,t w ) = 2c- 2t l Q (t-t w )[l + l 1 ]{t + t w )-e- 2it+t ^[l Q + li] 2 (t + t w ), (43) 
Xd(Mw) = 2e- 2t {[l Q + l 1 }(2t)-l Q (t-t w )[I + l 1 }{t + t w )}. (44) 

These results can be written in a more physically intuitive way in terms of the concentration of domain 
walls c(t), Eq. (|3 II) . and the return probability p r (r) = e _T Io(r) of a continuous-time random walker 
on a discrete, one-dimensional lattice 0|- Expressing all time dependencies in l|43(l . I|44(l via c(t) and 
Pi(t) yields the exact identities 

C d (t,t w ) = 4c(^-)[p r (t-t w )-c(^)], (45) 
Xd(Mw) = 4[c(t)-p t (t-t w )c(*&*-)]. (46) 

The fact that we find random walk-related quantities does not come as a surprise given that there 
is an exact mapping of zero temperature Glauber dynamics in the Ising chain to a diffusion-limited 
pair annihilation (DLPA) process 0|. The mapping follows by assigning to each bond + 1) the 
'particle' occupation number = ^(1 — SjS,+i) G {0,1} which signals the presence or absence of 
a domain wall. Glauber dynamics for the spins corresponds to independent random walks for the 
particles and coalescence of domains of aligned spins yields particle pair-annihilation. 

It follows from the definition 1)14(1 of the defect autocorrelation that Cd(t,t w ) — 
4 ({bi(t)bi(t w )) — (bi(t))(bi(t w ))) in fact also describes the particle autocorrelation in the DLPA pro- 
cess. We note that Il45fl is a non-trivial result. Assuming as in Ref. (3(| that the autocorrelation of the 
fraction c(t) of particles that still exist at time t is given by p r (t — t w ) and that these particles are un- 
correlated with the fraction c(t w ) — c(t) of particles that have disappeared via annihilation, one would 
conclude (6i(t)6i(tw)) = c(t)p T (t — t w ) + c(t)(c(t w ) — c(t)) and hence Cd(Mw) = 4c(t)(p r (t — t w ) — c(t)). 
This obviously differs from the exact solution (|45|) . As an approximation it holds for At <§; t w , but 
breaks down for At 3> i w where (|45ll yields Cd(t,t w ) « 2t w /(irAt 2 ) whereas the approximation gives 
Cd(t,t w ) ~ (V2— l)/(irAt). This shows that two-time correlations in C d (t,t w ) build up via a rather 
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subtle mechanism, the explanation of which in terms of DLPA would probably require knowledge of 
the inter-particle (i.e. domain size) distribution. Similarly, it appears that the result l|46|l cannot be 
obtained in a straightforward way. 
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FIG. 3: Defect autocorrelation (left) and susceptibility (right) versus At for waiting-times i„ 
10 1 , 10 2 , . . . 10 6 . Increasing waiting times correspond to decreasing values in the plot for small At. 



Now we turn to the dynamics of C d {t,t w ), Xd(t,t w ) - examples of which are shown in Fig. Bi- 
as given by l|43)l . i|4T)|) . The equal-time value of C d (t w ,t w ) — 4c(i w )(l — c(t w )) « 4c(i w ) « 2/ '\/irt w 
for t w 3> 1 decreases with t w , reflecting the decreasing number of particles in the DLPA process 
(or domain walls in the spin chain). In the regime At <C t w the two-time correlation C d (t,t w ) ~ 
4c(t w )p r (At) drops from its initial value due to the random walk motion of the particles around 
their initial positions at t w . In the aging limit of large At and t w one has the expansion C d (t, t w ) ~ 
2/(n^/t + tZ) (l/Vt^t^ - 1/Vt + tZ). This crosses over from C d (t, t w ) ~ 2/(ny/2At t w ) for At <C 
t w , where it connects smoothly to the initial drop for At of 0(1) since p r (At) w l/\/27r At for large 
At, to C d (t,t w ) ss 2t w /(TrAt 2 ) for At ^> t w . The integrated response Xd(Mw) is non-monotonic 
in At and increases on an 0(1) time scale in At from its initial value Xd(t w ,t w ) = to a plateau 
Xd(t, tw) ~ 2/yJnt w for At <C t w according to Xd(t : t w ) ~ 4c(t w )(l — p T (t — t w ))- This crossover is clear 
from the spin-chain dynamics: the perturbation associated with Xd(t,t w ) is 5H ~ —hsiSi+i which 
simply increases the coupling between sites i, i + 1. This enforces alignment of the spins Si and Si+i 
and hence increases (si(t)si+i(t)) on a microscopic time scale. In the aging limit the leading term in 
the integrated response is just Xd(t,t w ) ~ 2/y/iri which connects to the plateau Xd(t,t w ) « 2/y/nt w 



for At <C t w but eventually decreases as Xd{t,t w ) « 2/ynAt for At ^> t w . 

For constructing an FD plot (Fig. we are interested in keeping t fixed and varying t w between 
and t; the functions C d (t, t w ) and Xd(t>t w ) are then monotonic in t w . In fact the exact expressions 
(05), EH) satisfy C d (t, t) - C d (t, t w ) = Xd(t, i w ) + 4[c 2 ((i + t w )/2) - c 2 (t)}. Dividing this relation by 
the equal time value C d (t,t) yields the relevant normalized quantities 

1 - C d (t, t w ) = U t, M + ^fr^f - fcfc M + ^ j • (47) 

In the limit t — > oo the extra term in l|47(l vanishes and we get 1 — C d (t,t w ) — Xd(t,t w ) for all 
< t w /t < 1. This, however, does not imply that equilibrium FDT holds, i.e. X d (t,t w ) = 1. In fact, 
working out (d/dt w )C d (t,t w ) and R d (t,t w ) from 1)4*5) . j[""|") and expanding their ratio in the aging 
limit gives 

X d (t,t w ) M* + *w) _^ =. (48) 

iw(^ + tw) + (t-tw)v^ r ^ 

The FDR (|48|1 is a function of the time ratio t w /t and crosses over from X d (t, f w ) = 1 for i w /i — > 1 to 
Ad(t,i w ) = X^ — for t w /t — ► (Fig. 0). This seemingly paradoxical result can easily be explained 
in terms of the expansions given above. In the regime At <C t (which is equivalent to At <C t w , 
as considered before) we have, up to subleading corrections for t — > oo, C d (t,t w ) w 4c(i)p r (t — f w ) 
and Xd(t;^w) ~ 4c(t)(l — p r (t — t w )). So equilibrium FDT indeed holds in this regime and the 
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1-CAf ) 

FIG. 4: Normalized FD-plot (left) and the corresponding FDR versus tw/t (right) for the defect autocorre- 
lation and susceptibility. In both plots t is kept fixed, giving a one-to-one correspondence of the curves, and 
varies over the range t = 10° (dotted), 10 1/2 , 10 1 , 10 3/2 , 10 2 (solid). The curves of X d (Mw) for t = 10 3/2 , 10 2 
are almost indistinguishable and very close to the limit curve l@Ht - 



DLPA proces s is, to leading order, just an ensemble of independent random walks. Now recall that 
p T (At) « l/v27rAt for At ^> 1. So at the point where this approximation breaks down, At w t, the 
value of, e.g., Cd(t,t w ) decreases to an arbitrary small fraction of Cd(t,t) as t increases. This leads 
to a straight line segment which eventually covers the whole of the (normalized) FD plot while the 
size of the non-trivial region shrinks as 1/Vt. In the latter part, for Cd(t, t w ) <C one has from 

the aging expansions of C<j(i, i w ) and Xd(t, t w ) 

1-Xd(*,tw)w -^= + ^Cj(t,t w ). (49) 

Hence the FD plot indeed turns horizontal as Cd(t, t w ) approaches zero, consistent with 148(1 . In 
summary, an FD plot is not the appropriate representation for the FDT violation measured by the 
defect autocorrelation and response. A plot of the FDR as a function of t w /t, however, converges to 
the non-trivial limit curve given by ((48(1 as times diverge, see Fig. 0] 



2. Uniform field: Coherent functions 

For uniform covariances q n — 1 the defect observable Od is equivalent to the total energy of the 
system. According to 1(23(1 we have C e (t,t w ) = C(0;t,t w ) which may be simplified to give 

C„(t, t w ) = 4e- 2 *[I + Ix](2t) - e- 2 (* + *»>[3Io + 4I a + I 2 ](2(t + t w )). (50) 

This result again has an analog in the associated DLPA process, where it describes the normalized two- 
time correlation of the total number ofparticles A/", C e (t,t w ) = (4/A) [{Af{t)J\f(t w )) - (J\f(i)){J\f(t w ))]. 
A result similar to 1(50(1 was given in [44j , for initial conditions corresponding formally to equilibrium 
at inverse temperature 1/T = — oo. Up to a factor of 4 which appears to be missing in |44j . it coincides 
with (|5U|l for large t w , where one finds the simple scaling form C e (t,t w ) ~ 4/Y / 7r(l/v / t — 1/ \/i + t w ). 
At equal times, fluctuations in the energy follow as C c (t w , t w ) ~ (2 — v / 2)Cd(^w, i w )- This shows that 
in C c (t w ,t w ) = 5Z„ C„(t w ,t w ) the non-local (n ^ 0) terms make a contribution — (V2— l)Cd(t w , t w ), 
of the same order as the local term but with opposite sign. For At -C t w the two-time correlation 
C c (t, t w ) « C c (t w , t w ) has a plateau but it decreases as C e (t, t w ) ~ 2t w / '(Aty w At) when At 3> t w . 

By setting fc = in the Fourier transform l(B4(l we find that Xe(t, t w ) — x(0; t, t w ) = at all times. 
This is for the simple reason that the perturbation is proportional to the Hamiltonian and therefore 
just rescales the temperature, which obviously has no effect in the T — > limit considered here. We 
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note that Xe(t,t w ) — J2 n Xn{t,t w ) = implies that the sum over all cross-susceptibilities (n ^ 0) 
exactly balances the local susceptibility Xd(t, ^w) = Xo(*, *w)- 

An FD-plot for the energy is obviously just a horizontal line and the corresponding FDR is 
X e (t,tw) = X£° = = 0. This matches our findings in Sec. Ill Bl in the sense that the FD-plot for 
the coherent observable is a straight line whose slope is the X°° of the incoherent observable. 



3. Short-range correlated field 

We have seen above that for local defect observables the FD plot is not appropriate for determining 
FDT violation effects, since it converges to a straight line in the aging limit, ft turns out that the 
same holds for defect observables defined by short-range correlated fields. To see this, recall from 1)13(1 
that e.g. the correlation function C(t, t w ) of the observable is a weighted sum of the non-local defect 
correlations, and focus on the regime At = 0(1) that dominates the FD plot for large £ w or t. 
From 1(20(1 . one then easily shows that whenever a non-local term with given n ^ is of the same 
order as the local contribution Cd(t,t w ) = Co(t 7 t w ) ~ 2/-y/7rf w e~ A *Io(At), it can be written as 

C n (t, t w ) ~ -^=e- A *I„(At). (51) 

In the same regime the expression for the non-local susceptibility Xn(t,t w ) is identical apart from a 
minus sign. This shows that, whatever the short-ranged field correlations q n , the FD plot of x{t>tw) 
vs C(t, t w ) for the observable considered becomes trivial for long times, just as in the case q n — 8 n> o. 
We therefore focus on the FDR in the following, which requires analysis of (d/dt w )C(t,t w ) and the 
response R(t,t w ) = — (d/dt w )x(t, i w ) an d should become non-trival in the aging limit. 

By analogy with the results presented in Sec. Ill B 31 for spin observables, we will show that the 
FDR becomes identical to that for the incoherent functions in the aging limit. The procedure is 
again to prove that the Fourier transforms of the defect functions (d/dt w )C(k;t,t w ) and R(k;t,t w ) 
are representations of S in the aging-limit and with appropriate normalization. 

The expressions l(B3(l for C(k;t,t w ) and the one that follows from l(B4(l for R(k;t,t w ) are rather 
complicated and it is a priori not clear how they behave as times diverge. Asymptotic expansions in 
the aging limit t, t w — > oo with e < t w /t < 1 — 8 fixed (e, 8 > 0) and \k\ < K where K = c/y/t^ (c > 
arbitrarily large but finite), however, capture the relevant features of C(k;t,t w ), R(k;t,t w ) and have 
a considerably simpler form: 

C(k;t,t w ) ~ -ifI e - fc2 (^)/W-^e- t2 («-)/ 4 } (52) 
V 71 " L v t v* + *w J 

+ 2fce- fc2 ( t+t »)/ 2 {Erfi (fc^f ) - Erfi (fc^±S) J , 
R{k - 1 , tw ) ~ -L ( hL) k 2 e - fc2 (' 2 -*w)/(«). (53) 



Erfi(x) is the error function with imaginary argument: Ern(x) = (l/i)Erf (ix) |40|. Note that the 
arguments of all exponentials and the Erfi's are of 0(1) if t w /t and k are in the specified range. For 
|fc| larger than 0(l/y/t^) the results (|52*|l . f53^l do not apply. 

In l(52|) the growth of Erfi(x) ~ e x /(t/ttx) is over-compensated by the exponential prefactor and 
so we can make C(K;t,t w ) arbitrarily small by increasing c. For larger k, \k\ > K, the values of 
C(fc; t, t w ) as given by l(B3(l also turn out to be insignificant. Therefore C(k\ t, t w ) develops an infinitely 
sharp peak of width 0(l/y/t^) at k = in the aging limit and becomes a realization of S(k) when 
normalized by Cd(i, tw), in analogy with ((33(1 . Differentiating 1(52(1 w.r.t. i w turns out to reproduce the 
rigorous expansion for (d/dt Vf )C(k;t,t vr ) and similar arguments apply. Hence C(t,t w ) ~ N c Cd(t,t w ) 
and (d/dt w )C(t,t w ) ~ N c (d/dt V! )Cd(t,t w ) for any short-range correlated field with N c ^ 0. 

The expansion ((53(1 for the response function R(k;t,t w ) also peaks sharply in the region |fc| < K 
near k = 0; it follows from l(B4(l that outside this fc-range R(k;t,t w ) is insignificantly small again. 
We have, however, R(0;t,t w ) — at all times. Nevertheless, Rd(t,t w ) yields normalization and 
the ratio of both vanishes in the aging-limit for any k ^ (modulo 27r). So R(k;t,t w ) becomes 
a realization of S(k) when normalized by Rd(t,t w ), i.e. two infinitely sharp peaks at + and - . 
Therefore R(t,t w ) ~ N c R^(t,t w ) for any short-range correlated field with iV c ^ 0. 
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Since R(t,t w ) ~ N c Rd(t,t w ) and (d/dt w )C(t 7 t w ) ~ N c (d/dt w )Cd{t,t w ), any defect observable Od 
with short-range correlated fields ti and A c 7^ ultimately gives the same FDR as the incoherent 
functions. The scaling of the peaks in l|52|) . I|53[l implies associated time-dependent length scales in 
real space. As in the spin observable case the FDR will thus display a crossover (see Fig. 01 when 
these length scales become comparable with the length over which the fields e; are correlated. Wc 
note finally that, in contrast to the response R n {t,t w ) discussed above, the integrated response or 
susceptibility Xn(t, t w ) displays somewhat unusual behaviour; e.g. the local value Xo(A *w) dominates 
the non-local terms for all times, so that to leading order there is no real-space length scale associated 
with the defect susceptibility. One also finds non-trivial features in the FD plots and FDRs for the 
cross-correlations C n (t,t w ) and susceptibilities Xn(t,t w ) |if|. However, in the aging-limit and for any 
fixed n^O the FDR for the local observable (n — 0) is recovered as in the spin observable case. 



4- Infinite-range correlated field 



We next consider the FDR for observables defined by infinite-range correlated fields. As for short- 
range correlated fields, we will not discuss the integrated quantities C(t,t w ) and t w ) in detail. 
One finds again that these give a trivial FD plot for long times, although the argument for this is 
somewhat more subtle than in (|51|l because one needs to consider an infinite range of distances n. 

The exact expressions for the two-time correlation functions and susceptibilities for defect ob- 
servables with power-law covariances follow from l|23[) by substitution of C(k;t,t w ), Eq. (|B3|) . 
x(fc;i, t w ), Eq. ljB4|) . and qp{k), Eq. The response R(t,t w ) is then obtained from x(Mw) by 

R(t,t w ) = ~(d/dt w )x(t,t w ) as usual. The resulting equations are rather bulky and too complex for 
a meaningful discussion. So we immediately turn to the aging limit, where we can use the following, 
asymptotically exact, approximations. Firstly we replace the exact expressions in (I23[) for C(fc; t, t w ), 
R{k;t,t w ) by $5$. Although do not hold outside the range k G [-K,+K], the con- 

tributions to the fc-integrals are subleading. Secondly, as the integrands have infinitely sharp peaks 
at k = in the aging limit, we may replace qp(k) by the leading term of its expansion at k = 0, i.e. 
replace sin(fc/2) by k/2 in H25|) . This in turn allows us to extend the limits of integration in 
from — 7r, +7r to —00, +00, whereby we again just accumulate subleading errors. Having made these 
approximations, which still yield asymptotically exact results, the k integrations can be evaluated 
and we get 



Cp(Mw) ~ ~r(^){*^ (f -tiy 



where we have introduced the short-hand 



F(a;^)=l-(l-a)2- i ^(B(I ) l-f;^)-B(I,l-|;|)) 



(54) 
(55) 

(56) 



In the limit a — * 1, Eqs. (|54|l . I|55(l reduce to the asymptotic expansions of the incoherent functions 
Cd(t,t w ), Rd(t,t w ) while a — ► gives the asymptotic expansions of the coherent ones, i.e. C e (t,t w ) 
from and R e (t, t w ) = from JSHJ). So the power-law covariances lf2l))l again allow us to interpolate 
between local and global observables. For intermediate exponents < a < 1 the two-time correlations 
in Od decrease as t^, 1 ^ 2 At~ a / 2 in the regime 1 < A( <C t w and cross over to t w At~( 3+a ^ 2 for 
K t w < At. The response R P (t, t w ) behaves as Q 112 At-^+^Z 2 for 1«A(«( W and t w At- ( - 5+a ^ 2 
for 1 t w At. An aging expansion for the FDR again gives non-trivial curves. The derivative 
(d/dt w )Cp(t,t w ) follows correctly by differentiating the expansion (|54|l which, together with l|55|) 
yields 



Xp(t,t v 



1 



1 



2o 



2a 



t 



t + t v 



t-t v 



(57) 



The FDR Xp(t,t w ) is a function of t w /t only and interpolates between the FDR (|48|l for the local 
defect observables (a — > 1) and X c (t, t w ) — for the energy (a — > 0). For any power < a < 1 l|57|) 
crosses over from Xp(t,t w ) = 1 for At <§; t w to Xp(t } t w ) = = for At ^> t w (see Fig. EJ). 
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FIG. 5: Left: Time evolution of the defect-observable FDR for a random field with Lorentzian covariances 
12411 and a = 10 2 , correlated over N c = a7rcotha7r « air sites. For each curve t„ is kept fixed, varying over 
the range i w = 10°, . . . 10 6 , oo from bottom to top. The curves for t w = 10°, 10 1 , 10 2 are flat; in this regime 
the observable is effectively identical to the energy. For t w = 10 3 , 10 4 , 10 s , 10 6 we see the crossover to the 
limit curve for t w — > oo given by 148H and corresponding to the incoherent observable. Right: Limit curves 
of the FDR versus t w /t for t, t w — » oo. From bottom to top these correspond to power-law covariances with 
exponent a = 0.1,0.2,0.4,0.6,0.8. The top curve shows the random-field limit a — > 1 for power-laws, and 
also applies to any short-range correlated field Sec. Ill C 31 with iV c / or, in the extreme case, the incoherent 
functions of Sec. Ill CT1 

5. Harmonically correlated fields and X°° 

In contrast to spin observables it appears that for defect observables we generically find X°° = 0. 
To prove this claim we may again follow the approach presented in Sec. lIIB5l Introducing an FDR for 
defect Fourier modes X(k:t,t w ) according to (|40|l based on the two-time defect correlation function 
(|B3p and susceptibility (|B4|I allows us to write the FDR for any defect observable in the form (|4ip. 
The full expression for X(k;t, t w ) is rather complicated and, in contrast to l|42|) . retains a non- 
trivial dependence on k, t and t w . The only general features are X(0;t,t w ) = 0, since X(0-t,t w ) = 
X e (t,t w ) = 0, and X(±Tr;t,t w ) = 1 + 0{\ft^e~ ii "~) being independent of t and close to one for 
t w ^> 1. For intermediate values < \k\ < ir the FDR X(k;t,t w ) can in fact take arbitrarily large 
values for appropriate t w and t. To repeat the argument of Sec. Ill B 51 however, we just have to 
be able to take the limit t — > oo for fixed and finite t w . A scaling analysis of (d/dt w )C(k;t,t w ) 
as obtained from l|B3(l shows that this quantity develops an infinitely sharp peak at k = 0. Hence 
the normalized distribution of q(k)(d/dt w )C(k;t,t w ) over — 7r < k < tt becomes, for sufficiently well 
behaved functions q(k), a realization of 6(h) as t — > oo. Remarkably, the FDR for defects X(k; t, t w ) 
approaches a simple, smooth function on tt < k < ix in the same limit 

lim X(k; t, t w ) = 2 sin 2 \ ( 1 + 4t w cos 2 ^] (2 - e" 4 *" si " 2 i) . (58) 
t^oo 4 y 4 I V / 

Together, these two facts imply that out of the spectrum of FDRs X(k;t,t w ) for Fourier modes 
k, the long time limit < -> 00 again selects the contributions associated with infinite length scales 
(k — 0). These are, in the limit, given by l|58(l and equal to zero. The FDR (|41|l for defect observables 
X(t, i w ) — > thus vanishes as t — > 00 regardless of the choice of q(k), except in pathological cases as 
discussed in Sec. Ill B"5l Note that because (|58f) for k = gives a vanishing result for any t w , one in 
fact has lim^oo X(t, t w ) = X°° — 0, without needing to take t w — > 00. 

We note finally that the behaviour at short wavelengths is rather more complex for defect observ- 
ables than for spins. In particular, even for what one might expect to be 'equilibrated' wavelengths, 
k 2 t w ^S> 1, it is not true that X(k;t,t w ) w 1 for all times t, and X deviates significantly from this 
simple value for large time differences At ^> t w as can be seen from (|55)l . 



18 



D. Physical discussion 



We saw above that apart from pathological exceptions all spin and defect observables give identical 
values for the asymptotic FDR X°° , with X°° = 1/2 for spin observables and X°° = for defect 
observables. These slopes are most easily read off from the FDT plots for the coherent observables 
(magnetization and energy, respectively) which become straight lines in the long-time limit. 

It is natural to ask how these results would extend to observables other than those we have consid- 
ered, such as O = €iSiSi+2 which involves spin pairs at distance two. We have worked out explicitly 
the FD properties based on the general solutions given in |3Sj for the coherent and incoherent versions 
for this observable |45| ; one finds that they are, up to subdominant corrections, identical to those for 
O = 2 ejSjSj+i. The physical interpretation is simple: = —1 if there is exactly one domain 

wall between spins i and i + 2, while SiSi+2 = 1 if there is no domain wall or if there are two. The last 
alternative, however, is suppressed in the aging limit where typical distances between domain walls 
scale as and so s^Si+2 ~ SjS-t-i-i + Sj+iSj+2 — 1- For the coherent observable (q n = [eje$+ n ] = 1), 
this directly explains our observation; for the incoherent version (q n = 6„ t o) it follows from the fact 
that the correlations of s^Si+i and Si + iSi+2 are identical to the autocorrelations of SiSi + i in the aging 
limit. 

By a similar reasoning we can now predict the FD behaviour of higher-order observables of the 
form 



with k > 3 and j\ = 0; j2, ■ ■ • ,jk specify the relative displacements of the spins in the fc-th order 
products. For even k, each term Si+j,, again has a sign depending on the number of domain walls 
between spins i and i+jk- In the aging limit, configurations where more than one domain wall occurs 
can be neglected, so that we can replace the product by SiSi + j k . By the same argument as above, 
this is essentially equivalent to jk X SiSi+i and so should again give X°° = 0. 

For odd k, on the other hand, the sign of JT Si + j v is essentially determined by the sign of the 
domain which the spin Sj finds itself in. The leading contribution is now given by configurations 
with no domain walls between Sj and Si+j k . Configurations with at least one domain wall are again 
suppressed in the aging limit. We can therefore replace the product simply by Sj to leading order, 
giving an asymptotic FDR of X°° — 1/2 as for genuine first order spin observables. 

The fact that observables of even and odd order behave in different ways can also be motivated 
mathematically from the hierarchy of the equations obeyed by the multi-spin correlation functions 
where the even and odd orders k turn out to decouple completely (ssj. This is a peculiarity of the 
one-dimensional Ising model, whereas in the generic case one would expect all levels of the hierarchy 
to couple to each other, resulting in a unique value of X°° . We indeed find strong evidence for this 
in the two-dimensional case below. 

From a more physical point of view, the existence of two different values of X°° could be related 
to the fact that in the one-dimensional chain at T = one has both a critical point and an ordered 
phase. The result X°° = for defect observables could thus be related to the ordinary results for 
coarsening in d > 2 after a quench to an ordered phase, while X 00 — 1/2 for the spins would reflect 
the critical aspects of coarsening at T = T c = 0. 

It mi ght be interesting - though rather complicated - to use the methods described above and 
in Ref. [33] to study higher order observables different from (|59|) for which the above leading order 
approximations do not apply, for example, 1 — SiSi+i — Si + is.; + 2 + SiSj+2, which corresponds to a 
quadratic operator in bond variables, 46j&j+i. We are currently exploring this issue. 

To recap, the central result of this section is that (almost) all observables of the form l|59[) interpolate 
between an equilibrium like behaviour with X = 1 and an asymptotic FDR X°° . The latter are given 
by the values of X(k —> 0), as was argued in j^J IS EH- We have shown that the most efficient 
way of extracting X°° is by studying coherent functions. These results motivate the following section 
where the 2c? Ising model is studied at criticality. 



k 




(59) 
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FIG. 6: FD plots for spin autocorrelation and response. Three waiting times, t w = 43, 179 and 460 are rep- 
resented by squares, circles and triangles, respectively. The dashed line with slope one shows the equilibrium 
FDT. The full and dash-dotted lines have slopes X£° = 0.34 and X£° = 0.26, respectively; these are discussed 
in the text. 

III. THE 2d ISING MODEL 

In this section, we report on numerical simulations of the 2d Ising model. It is defined by the 
Hamiltonian 

n = -^s iSj , (60) 

(id) 

where the Sj (i — 1, • • • , N) are N Ising spins located on the sites of a square lattice with periodic 
boundary conditions and linear size L; the sum is over nearest neighbor pairs. We perform Monte- 
Carlo simulations using a standard Metropolis algorithm where the spins are randomly updated. One 
Monte Carlo step represents N attempts to flip a spin. 

The system is prepared in a random state, corresponding to an infinite initial temperature. It is 
then quenched at t = to the critical temperature T c = ^ ln(l + y/2). As stated in Scc. lIDI we focus 
on the four natural FD relations for the Ising model, constructed from the coherent and incoherent 
dynamical functions of spin and defect observables. The system size we use is different for coherent 
and incoherent objects. Incoherent objects reflect the behaviour of individual spins or defects, and 
simulating a very large system is advantageous in that it makes an average over many initial conditions 
unnecessary. Coherent objects like C c (t,t w ) or C m (t,t w ), on the other hand, have an amplitude of 
order l/JV. One should thus simulate many initial conditions of the smallest possible system, with the 
opposite constraint that the system has to be out of equilibrium even for the largest simulated time 
scale, tsi mu , giving the condition £(t s imu) -C L. Our results are obtained with i s i mu = 10 5 , L = 300 
for coherent functions, and L — 500 for incoherent ones. Only a few samples over initial conditions 
are necessary for incoherent correlation functions, while 1000 initial conditions were sampled for 
coherent ones. This is also the number of realizations necessary to get the four susceptibilities we 
have computed. 

We now describe our results, starting with spin observables and then turn to defect observables. 

A. Spin observables 

The two-time scaling behavior of the incoherent spin functions C s (t,t w ) and Xs{t, t w ) has been the 
subject of a number of publications, as described in Sec. II CI We refer to the references cited there and 
directly present in Fig. the FD plot for the spin autocorrelation and susceptibility. A very similar 
FD plot has been reported in Ref. [2j|, although a somewhat different susceptibility J Q W dr R s (t, r) 
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FIG. 7: Tentative extrapolation of the infinite time slope of the FD plots of Fig. The lines are only 
suggestive, indicating that an asymptotic FDR of X°° = 0.34 is compatible with the data in the regime 
C s (t,t) — C s (t,i w ) ~ 1. The different symbols have the same meaning as in Fig. 



was plotted there, so that the FD plot looks reversed compared to Fig. Otherwise, we find the 
features anticipated in Sec. II CI The FD plot is characterized by an initial part which follows the 
equilibrium FDT, corresponding to short, equilibrated length scales. For larger time differences, the 
FD plot deviates from the FDT in a non-trivial manner due to the non-equilibrated fluctuations at 
small wavevectors. In the limit of large time differences the FD plot has a nonzero slope A s °°, in 
contrast to the zero slope obtained below the critical point 0, @ ■ These features make the FD plot 
rather similar to the one obtained in d = 1 , see Fig. [21 (left) . 

The infinite time value for the slope of the FD plot for the 2d Ising model was estimated in Ref. [2^ 
as X s °° = 0.26. We recognize from Fig. El that the crossover from X = 1 to X£° < 1 takes place 
over a very small range of the correlator, and that a precise determination of the infinite-time value 
of the FDR is difficult. A tentative numerical extrapolation is shown Fig. [7| where the quantity 
(Xs(t,t w ) — X s (oo, t w ))/C s (t,t w ) is plotted against C s (t,t) — C s (t,t w ); as the abscissa approaches 1 
(i.e. for large time differences), the ordinate should converge to X^°. The figure shows that the value 
X£° = 0.34 is compatible with the data, but even though we use larger waiting times than in Ref. |22| 
there is substantial scatter in the points. However, we have more precise estimates of X£° to guide 
us, as we now describe. 

The study of the model in d = 1 in Sec. [H] showed that the crossover from X = 1 to X^° for 
spin functions reflected the different dynamics of large and small wavevectors which for the d = 2 
case would be defined according to fc£(t w ) 1. The dynamical behaviour of the small wavevectors 
was governed by the asymptotic FDR X%°. This suggests that a much simpler measurement of X£° 
should be possible by focusing on the k — ► limit, i.e. by measuring the correlation and susceptibility 
of the magnetization density m(i). The resulting FD plot is reported in Fig. El As for the Id case, a 
very simple result is obtained, with the FD plot extremely well fitted by a simple straight line. The 
fit holds for several decades of time, t w < t < i s imui I0r each waiting time t w that we have considered, 
providing strong evidence that X m (t,t w ) — X^ at all times. Furthermore, the slopes of the three 
curves in Fig. IHlare very close to one another, and this allows us to report the value 

X™ = 0.340 ± 0.005. (61) 

This is the value we used to fit the data for the incoherent spin functions in Figs. El and demon- 
strating that the data are consistent with the equality X^° = X^ . This is somewhat different from 
the value reported in Ref. |29| , but we believe that our measurement from the magnetization is much 
more reliable than the extrapolation of the incoherent spin functions, as explained above. We note 
also that this value is in extremely good agreement with the two-loop expansion value reported in 
Ref. [33l ]. However, unlike the Id case, we do not have a simple physical argument to explain the 
actual numerical value. 
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FIG. 8: FD plots of correlation and susceptibility for the magnetization. The three curves are for waiting 
times t v = 46, 193 and 720 (bottom to top). The curves for t w = 193 and t w = 720 have been shifted 
vertically for clarity, since they would otherwise overlap with the curve for t„ = 46; the unshifted curves all 
pass through the origin as they should. The dashed line is the equilibrium FDT. The full lines have slope 
X™ = 0.34. 




FIG. 9: FD plots for defect autocorrelation and susceptibility. Data for two waiting times, t w = 43 (squares) 
and t w = 179 (circles), are presented. The second one has been vertically shifted for clarity. The full lines 
represent the equilibrium FDT. 



B. Defect observables 



We now turn to defect observables. The simplest functions to consider are the defect autocorrelation 
function and the conjugate susceptibility. These quantities have been studied recently for kinctically 
constrained Ising models (in particular the Fredrickson- Andersen model in Id), where they were shown 
to give rise to simple FD plots |2^. We present the corresponding FD plot for the 2c? Ising model 
in Fig. El Again an apparently very simple result is obtained, with the FD plot very well fitted by 
the equilibrium straight line with X = 1. This is an unexpected result, since the system is far from 
equilibrium as was demonstrated b y th e study of spin observables in the previous section. It could 
also be taken to imply, as in Ref. [2f|, that the asymptotic value of the FDR associated with the 
defects has the equilibrium value X™ — 1. 
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FIG. 10: FD plots for the energy density. The three curves are for t w = 80, 193, 464 (bottom to top). The 
curves for t w = 193 and t w = 464 have been vertically shifted for clarity, and would otherwise again pass 
through the origin. The dashed curve is the equilibrium FDT; the full lines have slope = 0.34. 

Our above study of the Id model again clarifies the situation. There, we found that the inco- 
herent dynamical functions of the defects exhibited a crossover from equilibrium to non-equilibrium 
behaviour, but that the non-equilibrium part was barely visible in an FD plot since the crossover 
occurs when correlators have already decayed to very small values. This suggests that the apparent 
equilibrium behaviour observed in simulations for the 2d Ising and Id Fredrickson- Andersen models 
is simply a good approximation to numerical data, but may miss non-trivial FD relations at large 
times due to limitations in the numerical analysis. However, as for the spin observables, the solution 
to this problem is straightforward and consists in focusing on the k — > limit. We thus investigate 
next the coherent functions for the defects which are the autocorrelation and susceptibility for the 
energy density. 

The resulting FD plot for the energy density is shown in Fig. 1101 As for the magnetization, very 
good fits by pure straight lines are obtained, implying the equality X e (t,t w ) = X£° . Note, however, 
that these plots have more noise than the ones for the magnetization. This is due to the fact that 
the abscissa now involves a genuine connected correlator, in which the nonzero average of the energy 
density needs to be subtracted off. Nonetheless, the slopes of the FD plots in Fig. ^] are very close 
to one another and give the result 

X c °° = 0.33 ± 0.02. (62) 

An important outcome of this paper is that this value is compatible, within error bars, with the value 
reported above for the asymptotic FDT for the magnetization and the spins. This strongly suggests 
that the various infinite time FDRs that we have measured in the 2d Ising model are all equal, and 
that the critical point of the 2d Ising model is described by a single new universal quantity 

X™ = X™ = X™ = X c °° eI=°« 0.340. (63) 

IV. SUMMARY AND DISCUSSION 

In this paper we have studied the relation between two-time, multi-spin, correlation and response 
functions in the non-equilibrium critical dynamics of Ising models, analytically in the d = 1 case, 
and numerically in d = 2. We have shown that FDRs, while observable-dependent, fall into well- 
defined classes, which are qualitatively similar to those observed in various glassy systems. All FDT 
violations can be understood by considering separately the contributions from large wavevectors, 
which are at quasi-equilibrium and obey FDT, and from small wavevectors where a generalized FDT 
holds with a non-trivial fluctuation-dissipation ratio X°° — X(k — > 0). In d = 1, we find through 
exact calculations X°° = | for spin observables and X°° = for defect observables. In d — 2 we find 



23 



numerically a unique X°° ~ 0.34 for all observables. These results suggest that the definition of an 
effective temperature T c s = T/X°° for large length scales is generically possible in non-equilibrium 
critical dynamics. 

Further, this work also suggests many interesting lines for future investigation. An important 
question is what are the limiting FDRs in diffusive models that are analogous to the Id Ising model 
but have glassy features, for example the one-spin facilitated Fredrickson-Andersen model |47j or 
symmetric plaquette models |3fij |. Also, it would be interesting to confirm our results for the 2d 
Ising model by analyzing higher order correlation functions by means of the RG techniques used in 
Ref. [33L I34I l35| to confirm the uniqueness of the FDR. This would make this function an interesting 
quantity to study in more generic non-equilibrium situations such as driven interfaces or driven 
diffusive systems. 
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APPENDIX A: MODIFIED BESSEL FUNCTIONS 



Here we briefly summarize the main properties of modified Bessel functions l n (x) that are relevant 
for the analysis given above. A comprehensive description may be found in |40| . For integer order n, 
I„ (x) has the integral representation 



TT 

l n {x) = / — cos nip e 
J n 



(Al) 



from which the functional relations 



d 1 2n 

—I n (x) = -[I n -i+I n+1 ](x) and — I„(a;) = [I„_i - I„+i](a;) (A2) 

follow immediately. In particular it is clear from (|A1|) that I_ n (ir) = l n (x) and I n (— x) = (— I)"I„(a;). 
The aging expansions of our results are based on the asymptotic formula 

Ux) = ^(l+ 1 -^+o(±)), (A3) 



. 2irx V 8x 

which applies in the limit of large arguments x for fixed order n. For the derivation of the Fourier 
transforms of multi-spin correlation and response functions we use 

^V m;E i n (a) = e acosx , (A4) 

n 

k 



^e mk l n - m {a)l n+m (a) = I 2m 2acos-j, (A5) 



2 

^e- mfc I„(a)(I n+m (fe) + ln- m {b)) = 2T m { b + a ™ sk \ Im(A)) (A6) 

n ^ ^ 

where l|A5|l . (|A6|I follow from (|A1|I . the well known identity l|A4(l and trigonometric relations. In 



[6j < a < b is required, A = yja 2 + b 2 + 2ab cos k and the T„(x) = cos(narccosx) are Chebyshev 
polynomials of degree n in x. 
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APPENDIX B: FOURIER TRANSFORMS 



The Fourier transforms of spin correlation and response functions l)18|). (|19|l follow immediately 
when using ilA4(l : 



C(fc;t,tw) = e -( t+t -)( 1 - cosfc ^ 1+ / dr e - rcosfc [I + I x ](t) I , 



X(A;t,*w) = IJdr e-^ 1 -^ e- 2T [l + 2l 1+ l 2 ](2T). 



(Bl) 



(B2) 



For defect correlations, however, a direct transformation of (121)11 yields a rather intractable expression. 
Therefore we first rewrite H2()(l using the identity (which can be verified by differentiation) 

t + tw 

e-(*-*-)l„(t - t w ) = e-(* +t -)l n (i + tw) - /" dr -e^I^x - 2I„ + I„+i](r), 



as 



t+t„ 



C n (t,t w ) = e- 2(t + t ^[ll-l n ^ 1 l n+1 }(t + t w ) + ^ J dre-^+^x 

t— t w 

{[I„_i - I„+i](r) [I„_i - I„+i](t + tw) - Pn-i - 2I„ + I„+i](r) [I n _i + 2I„ + I n+1 ] (t + t w ) } • 

Now, utilizing (|A2(I and expressing factors of n as derivatives w.r.t. k, the Fourier series for C(k; t, t w ) 
may be written in the form: 



f 4 a 2 



t + tw 

J] ^ K -In-lln+l] (* + *w)"^ y d 



re (Ww+^x 



-(t + t w )<% 2 



^e- mfe I n (r)I n (t + t w ) + 2 ( — - 1 ) e- infe I»(r)[I„_i + 2I n + I„ +1 ](t + t w ) 



All summations in this expression can be evaluated via l)A5|l . (|A6|1 . Some fairly complicated algebra 
is required to simplify the resulting expression, but finally one obtains the compact result 



C{k\ t, t v 



- 2 (t+t w ) Ii(2(t + t w )cos|) 



(f + t w ) COS ■ 



(B3) 



4 / dre" 2 ^^ 



ll 1 (2A) + 2^sin 2 



I 1 (2A) + t -^h(2A) 



where A = y (t + t) 2 cos 2 (fc/2) + (t w — t) 2 sin 2 (fc/2). Eq. (|B3|) is the most convenient representation 
for C(k;t,t w ), both for numerical and analytical purposes. The calculation of the Fourier transform 
of the defect susceptibility (|2"TJ) is comparatively easy; from (|A5j) . (|A6|) one finds 



X(fc;Mw) =2e" 2 N [I + Ii](2t)-I 



(2A) - ^ s 2 (fc/2) + t w sin 2 (fc/2) ii(2A) | 



.4 



(B4) 



with A = y^ 2 cos 2 (fc/2) + tl sin 2 (fc/2). 



25 



APPENDIX C: POWER-LAW COVARIANCES 



Here we show that the covariances qp. n given in (|25|l follow a power-law as \n\ — ► oo and establish 
the link qp :ll = T~ 1 {qp(k)}. Let us first focus on the Fourier integral (|22H which - since qp(k) is even 
in k — may be written as 

r 2 (^) Jdk ( . ky- 1 

9P,n = \ / ^ sin- cos(nfc), (CI) 



r(a) J 2vr V 2 
o 

where < a < 1 as before. The simple substitution x = k/2 yields the solvable integral |4fj 

a x (—11" 

(sina;) Q ^ 1 cos2nx= : — ttj rr r . (C2) 



7T 



Using the functional relation [4£j B(x,y) = T(x)T(y)/T(x + y) for the Beta function ~B(x,y) and 
simplifying the remaining expression yields the result for qp, n given in (|25|l . Now we turn to the 
asymptotic behaviour of qp, n as \n\ — > oo. For n > 1 we may rewrite qp t m using r(x)r(l — x) = 
7r/ sin7ra; and T{x + 1) = xT(x), in the form 



f- 1 1 + a + 2k 

k—0 

It is obvious from (|C3|I that qp, n is monotonically decreasing and vanishes for n — ► oo as long as 
a > 0. It is equally clear that qp jn — 1 for a — > and <7p.„ = 5 ni o as a — > 1 (since gp !rl is even in n 
and qpfl — lVa). In order to understand the asymptotic behaviour of qp, n we take the logarithm of 
J3I) and use the bounds 

" n-l 

/ dfca fe < ^ a fc < a + / dfca fc , (C4) 
o fc=0 o 

which hold for any non-increasing function . For the case at hand the integrals can be solved easily. 
Exponentiating the result, multiplying by n a and taking the limit n — + oo then gives 



^\ TT ^ < < p.)-^^, (C5, 

which implies that there exists a finite constant c such that (7p.„ ~ cn~ a for < a < 1. 
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